Spin-Hall Conductivity and Electric Polarization in Metallic Thin Films 
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We predict theoretically that, when a normal metallic thin film (without bulk spin-orbit coupling, 
such as Cu or Al) is sandwiched by two insulators, two prominent effects arise due to the interfacial 
spin-orbit coupling: a giant spin-Hall conductivity due to the surface scattering and a transverse 
electric polarization due to the spin-dependent phase shift in the spinor wave functions. 
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Spin-orbit interaction, transferring angular momen- 
tum between electronic spins and orbital motion, has ex- 
tended the boundary of the field of spintronics towards 
a full-electric manipulation of spins without using mag- 
nets. By coupling the charge and spin currents, spin- 
orbit interaction has left its signature in bulk metalsi"— 
and semiconductorSii"— by the so-called spin-Hall effect;^ 
which catches much attention in academia and industry 
due to its interesting physics and potential applications. 

In noble metals such as Pt and Au, a strong bulk 
spin-orbit coupling drives the (inverse) spin-Hall effect, 
which has been thoroughly studied in numerous set- 
tings, providing an alternative to generate or detect spin 
currentFiSrJ^ Theoretical investigations, ab initio calcu- 
lations in particular, have strengthened our understand- 
ing of the role played by bulk impuritie o^^i^^ and doped 
surfaceSfi^ With a series of efforts to dope Cu using spin- 
orbit scatterers (see ReffT?! and references therein), a 
large spin-Hall angle has been reported in Cu thin film 
with bismuth impurities recentlyJ^ In this Rapid Com- 
munication, we propose an alternative method to realize 
the spin-Hall effect in Cu without doping but by sand- 
wiching the Cu thin film with two dissimilar insulators 
such as oxides or even vacua. The inversion symmetry 
breaking across the interfaces provides interfacial spin- 
orbit couplings (ISOCs), thus allowing metals such as Al 
or Cu to accommodate a spin-Hall conductivity that may 
be even larger than that caused by a bulk spin-orbit inter- 
action in noble metals. Meanwhile, we demonstrate that 
the ISOC also induces an in-plane wave- vector-dependent 
spatial separation of wave functions between opposite 
spin types along the confinement direction. Combined 
with the spin imbalance (between the majority and mi- 
nority bands) due to the structural asymmetry, such a 
spatial separation induces a transverse electric polariza- 
tion that is quadratic in the in-plane electric field. 

A schematic picture of the setup is shown in Fig. [T] 
The normal metal film has a constant thickness d along 
the transverse direction (z), while two interfaces are in 
the xy plane located at z± = zLd/2. The confinement of 



the insulators is described by finite potential steps, 

Vc - V+0iz - z+) + V-e{z- - z), 



(1) 



where V± is the height of the potential barrier at z± ~ 
zLd/2 and 6{z) is the Heaviside step function. Wc are 
interested in the surface scattering by the Rashba-type 
spin-orbit interaction^^ generated at interfaces between 
a normal metal film and insulating materials, such as 
metal oxides. The potential barriers Eq. ([T]) generate 
two electric fields Ec = — VVc that are localized at the 
interfaces and aligned oppositely to each other along the 
z direction, giving rise to a Rashba-type ISOC, 

„ X-V^d{z - z^} - X+V+d{z - z+) . „ 

LiR = [CTxPy - CTyPx) , 

(2) 

where X± is the spin-orbit coupling parameter for the 
corresponding interface, cr is the Pauli matrix, and p is 
the canonical momentum. 

Band structure and wave function. The full Hamilto- 
nian for this system is H = p^/ {2m) + Vc + Hb,. We first 
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FIG. 1. (Color online) Left: A metallic film of thickness d 
is sandwiched by two dissimilar insulating materials, such as 
oxides or vacua. The electric fields (solid blue arrows) at the 
interfaces point inward (±z), and the resulting Rashba mag- 
netic fields (dashed red arrows) point sideward (±t/) for an 
electron moving along the x direction. Right: The confining 
potential Vc in Eq. ([1} along the confinement direction (i.e., z 
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FIG. 2. (Color online) Left: Schematic view of the band 
structure in metallic thin film with interfacial Rashba spin- 
orbit coupling, majority (solid red) and minority (dashed 
blue). Right: Schematic picture of the spatial part of the 
wave-function for the first two subbands; the stepwise curve 
is the confinement potential profile and the three different 
curves are the wave function from the lowest-order perturba- 
tion Eq. (I3b|l (solid black for both majority and minority) and 
the exact wave function from Eq. 1)10^ (solid red for majority 
and dashed blue for minority) which shows the spatial sepa- 
ration of opposite spins. The discontinuity at the interfaces 
is due to the J-function-like potential in Hr. 



treat the ISOC with degenerate perturbation. Using 
the wave functions of a finite potential well and a stan- 
dard degenerate perturbation technique)^! assuming the 
energy splitting due to the ISOC is much smaller than 
the interchannel energy spacing, we obtain the energy 
eigenvalue and wave function for an eigenspinor labeled 
by spin polarization s = ±, an in-plane wave vector fc, 
and a transverse channel index n. 
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s—k\Xj^ — A_|_Eon , 
dp 



(3a) 
(3b) 



where tan6' = ky/k^, Eo=h'^TT'^ / {2mdl), and A is the 
area of the metal film. A schematic view of the band en- 
ergy and the spatial part of the wave function is shown 
in Fig. [3) The spatial part of the wave function ipnks 
in Eq. (j3bp is the same for majority (s = — ) and mi- 
nority (s = -f ) spins in the first-order perturbation cal- 
culation. The phase shift for both spins is given by 
S = mr/2 + {d+ - d^)mr/{2d) with d± = h/y/2niV± the 
approximate penetration depth into the barrier. By as- 
suming high confinement barriers (V±^Ep, given Ep the 
Fermi energy), the effective thickness de = d + dj^ + d-. 
Therefore, a finite potential well of thickness d can be 
viewed equivalcntly as an infinite well with thickness de- 
We shall point out that, although the eigenenergies for 
majority and minority are different, such a difference 
does not lead to any magnetism because of the exact 
cancellation among spins at different k directions. 

Spin-Hall conductivity. When A_|_7^A_, for each con- 
ducting channel (or transverse mode), the electron is 
characterized by a free Hamiltonian Hq^^ = /{2m) + 



Ei^n? augmented by an effective momentum-dependent 
Rashba magnetic field B^{p) = {2Xf /h){zxp) with 
a channel-dependent coupling constant X'^ = 2(A+ — 
X-)EQn'^ /de- For an electron in band n carrying an in- 
plane momentum p, assuming A+ > A_, the majority 
(minority) spin points to the {—B^) direction. 

When an electric field is applied along the x axis, 
the Rashba field B^'{p{t)) becomes time dependent 
since Px{t) = —eEx^ In the rotating frame that fol- 
lows B^'{p{t)), the time dependence is translated into 
a gauge field B^ — —ehpyEx/p^zi^ In the adiabatic 
limit {eEx'^Xj^ kp) , the majority (minority) spins align 
(antialign) with the total field JBj — B^ + B^, resulting 
in an out-of-plane (z) spin component 



(4) 



Consequently, in a metal film of volume V = dA, at the 
Fermi level, the spin current polarized in the z direction 
while flowing in the y direction is 



y yZ^ 



h Py eric ^ 



(5) 



Therefore, the spin-Hall conductivity caused by ISOC is 

'^kp. (6) 



e Uc 
8^~d 



where Uc — [kpd/irl is the total number of transverse 
channels. The symbol [aj denotes the largest integer 
that is smaller than a. Note that Eq. ([6]) is obtained for 
A+ > A_, the sign of cr^^ should change when A+ < A_, 
and cr^H = for A+ = A_. 

To extend the semiclassical picture above, we employ 
the Kubo formula to calculate the dc spin-Hall conduc- 
tivity (i.e., w — >■ 0) in the linear response regimes^ An 
electric current flowing along the x direction gives rise 
to a spin current that is polarized to the z axis while 
transporting along the y direction. 
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{Enks — En'ks')iEnks — En'ks' — fuj — irj)' 

(7) 

where fnks is the electron spin occupation number. Using 
the Heisenberg equation, the in-plane velocity operator 
is given by 



Jx{y) 



_ fi ^ X-V-S{z - z^) ^ X+V+S{z - z+) ^ 

- Ji '^■^W' 

(8) 

where the second term is the anomalous velocity due to 
the spin-orbit coupling that is localized at the interface. 
The deflnition of spin current polarized along the z di- 
rection is = (h/i){d-z,v}. Using the cigensolutions in 
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Eq. ([3]), WG have 

where ± = sgn(A+ — A_) and the last approximate value 
assumes a large number of transverse channels nc3>l. For 
two perfectly identical interfaces A_ = A+, the spin-Hall 
conductivity shall vanish. Equation ^ calculated from 
Kubo formula agrees with Eq. ([6]) derived semiclassically 
and comprises one of the main results of this communi- 
cation. 

To describe a spin-Hall system, one quantity often 
scrutinized in experiments is the spin-Hall angle defined 
as the ratio between cr^^ and the longitudinal conductiv- 
ity tJ^ , characterizing the efficiency of converting charge 
(spin) current into spin (charge) current in a (inverse) 
spin-Hall system. For a coherent ballistic conductor pre- 
sented in this communication, the dc longitudinal con- 
ductance (or resistance) measured in an experiment is 
dominated by scattering events at contactSF^^i Therefore, 
we expect the spin-Hall angle to depend on the specific 
geometry and material selection of the contacts. 

Transverse electric polarization. We now turn to an 
inviting effect, as caused by ISOC, that has not been dis- 
cussed before to our knowledge. The ISOC in Eq. ([2]) 
not only gives rise to an energy splitting for opposite 
spins, but also separates the wave functions in real space; 
i.e., majority and minority spins are shifted towards op- 
posite surfaces. This can be easily understood as the 
following: Assuming A+ > A_ , the potential barrier at 
surface z+ [z^) is decreased (increased) by a (5- function 
Rashba potential for majority spins; therefore majority 
spins tend to shift towards surface z+. On the opposite, 
the minority spins shift toward surface z_. Therefore, 
majority and minority spins are spatially separated along 
the transverse direction as shown schematically in Fig. [2l 
To quantify such an effect, the approximated wave func- 
tion in Eq. (j3bp is not enough, and we must seek an exact 
eigensolution to the full Hamiltonian H of the following 
form (inside the metal film): 

with normalization factor Cnks and the expression outside 
the metal is written similarly but with evanescent wave 
function in the z direction. In the limit of X±kp<^l and 

Qnks^-j- + [sk{X+ - A_) - {d+ + d_)] , (11a) 

T? TT HIT 

Snks^^ - ^ [sk{X+ + A_) - (d+ ~ d_)] . (lib) 

The spin-dependent transverse wave vector Qnks gives the 
energy splitting as in Eq. ([3a)) for the majority and mi- 
nority spins with Enks = /'^'^){lnks + ^^)- "^^^ spin- 
dependent phase shift S^ks means that the wave function 
for majority (minority) spins with s = — (-I-) is shifted 



toward negative (positive) z direction. The transverse 
shift is spin and k dependent, and can be quantified by 
the center of probability in the z direction: 

Zns{k) = j z\Tpkns{z)\'^dzKsk '^^ . (12) 

The last approximation omits a constant {d-\- — rf-)/2 
that does not contribute to the nonequilibrium prop- 
erties discussed below. For highly asymmetric inter- 
faces (i.e., A+^A_), the transverse shift \zn+{k) — 
Zn-{k)\(xrj{k)d, where rj{k) is the ratio between the in- 
terfacial Rashba energy splitting [third term in Eq. (|3a[) ] 
and transverse channel energy spacing [second term in 
Eq. ([3a|) ]. Therefore, by a realization of large ry, the ma- 
jority and minority spin channels can be spatially sepa- 
rated, and the spin-flip scattering is thus suppressed. 

We now discuss the effect of the transverse shift on 
the electric polarization. What interests us the most 
is the nonequilibrium response of the electric polariza- 
tion as a direct consequence of the ISOC. Because of 
the k dependence, the transverse shift z„s(fc) depends 
on the application of a current; therefore a nonequi- 
librium transverse electric polarization response to the 
in-plane electric field E can be calculated as Pz{E) — 
^Y.nksVris{k + 6k) - Zns{k)] with Sk = erE/h: 

The second expression is for a coherent conductor with 
an electron dwelling time r = l/vp and V = El. The 
quadratic dependence of £^ or is a result of symmetry 
in the in-plane dimensions. 

In an analogy to Hall effects, the electric polarization 
Eq. (fT3|) shall give rise to, across the confinement direc- 
tion, a voltage signal that is quadratic in the (in-plane) 
applied voltage. The electric polarization Eq. (fT3|l is 
a combination of two facts: (1) the spatial separation 
between majority and minority spins, which is propor- 
tional to A+ -|- A_, and (2) the spin imbalance between 
the majority and minority spins, which is proportional to 
A+ — A_. Equation psp comprises the other main result 
of this communication. 

We provide an estimate on the order of magnitude 
of the polarization density for the most asymmetric in- 
terfaces (e.g., A+^0 and A_ = 0). The existing value 
of Rashba parameter in literature ranges from = 
0.16 eVA2^ to an = 2.5 eV A,— depending on the mate- 
rials. Here, we take ur = 0.16 eV A, corresponding to a 
Rashba energy splitting A/? ~ 0.22 eV, for Cu with kp = 
1.36 A~^. We convert such an energy scale into A by 
AuRiAkpXEo/d. For a thin film of a thickness d = 10 nm 
and under an applied voltage, for example V = 100 nV, 
the polarization density Pz^2.1xl0~^'^ C m~^. Such an 
electric polarization is measurable (in a circuit enclosing 
two interfaces) under an ac in-plane electric field (or cur- 
rent): With a frequency / = 1 GHz. the magnitude of 
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the induced current density in the transverse direction 
is jz^dPz/dt^lO~^ A m~^. Meanwhile, the frequency 
of this induced current is doubled to 2/ because Pz de- 
pends quadratically on V. Such a frequency doubling, 
as a consequence of the symmetry of the present system, 
is qualitatively different from Pershin and Di Ventra in 
Ref. [2^, where the frequency doubling effect emerges 
from electron-electron interaction. In addition, we em- 
phasize the electric polarization in our study arises from 
the phase shift of a single-electron wave function along 
the out-of-plane z direction, while in Ref. [2^ the electric 
polarization arises from a nonlinear effect due to many- 
body electron-electron interaction and an inhomogeneous 
charge density. 

To the best of our knowledge, the transverse electric 
polarization induced by the ISOC predicted in this com- 
munication is a new effect that has not been discussed 
previously. This electric polarization manifests itself as 
a charge Hall effect and can be measured as a transverse 
voltage (or a current in a close circuit). Being quali- 
tatively different from conventional charge Hall effects 
driven by Lorentz force in normal metals and anoma- 
lous Hall effect in ferromagnetic media, the electric po- 
larization [Eq. (jl3p ] is quadratic (instead of linear) in the 
longitudinal electric field and does not require external 
magnetic field or ferromagnetism. 

Discussion and Conclusion. The spin-Hall conductiv- 
ity in Eq. ([9]) is derived for a ballistic sample where the 
bulk impurities are scarce. This is generally valid for thin 
film with thickness less than the mean free path and spin 
diffusion length. In an ultrathin film, surface roughness 
is a dominating scattering mechanismi^ In the present 
setting, the electric fields generated by the potential gra- 
dients are normal to interface while the presence of sur- 
face roughness effectively randomizes the field around an 
average direction that is still perpendicular to the av- 
erage interface. Therefore, we can summarize the total 
effect of surface roughness into an effective spin-orbit cou- 
pling Xi that is smaller than A^. Since the leading order 
spin-Hall conductivity docs not depend on A's, we argue 
that the spin-Hall effect survives the roughness scatter- 
ing. The seemingly abrupt jump in spin-Hall conductiv- 
ity [i.e., Eq. ([9])] at the point A+ = A_ is the manifesta- 
tion of neglecting electron momentum relaxation in such 



a ballistic conductor. In a realistic setup, the sandwich- 
type conductor is always connected to reservoirs where 
the electrons are relaxed and a smoother change shall 
appear. 

We need to point out that these results are, in many 
ways, different from the seminal works by Sinova et al. in 
a two-dimensional electron gas^ and Murakami et al. in a 
bulk semiconductor!^ First of all, in the present setup the 
spin-orbit coupling is neither intrinsic to the electrons (as 
in Ref. [ah nor arising from a particular band structure (as 
in Ref. |3), but due to the interface scattering. Second, 
our treatment to the spin-Hall conductivity is seemingly 
two-dimensional but the outcome highlights a bulk effect 
with a weak thickness (d) dependence, as long as the 
structure can be treated coherently. 

Third, and most importantly, because of the finite 
size in the z direction, such a sandwich-type struc- 
ture also accommodates a transverse electric polariza- 
tion (being qualitatively different from the spin polar- 
ization studied by Edelstein^i) through the ISOC, which 
is unique for the thin-film structure: The out-of-plane 
electric polarization along the z direction does not exist 
in two-dimensional systems considered, for example, in 
RefsHm and [13. 

In conclusion, we predict that in a coherent ballistic 
conductor that consists of an ultrathin normal metal 
(with negligible bulk spin-orbit coupling) film sand- 
wiched by dissimilar insulators, the Rashba-type ISOC 
supports a spin-Hall effect featured by a large spin-Hall 
conductivity that is independent of the interfacial Rashba 
coupling. The ISOC also causes a spatial separation in 
the transverse wave functions of different spin bands. In 
response to an in-plane current or electric field, such a 
spatial separation gives rise to a transverse electric po- 
larization that is quadratic in the in-plane field applied. 
The sandwich-type structure proposed in this communi- 
cation has potential applications to replace noble metals 
(such as Pt) as a source and detector for spin currents. 
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